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Abstract 

We study finite-time singularities in the linear advection-diffusion equation with a variable 
speed on a semi-infinite line. The variable speed is determined by an additional condition at 
the boundary, which models the dynamics of a contact line of a hydrodynamic flow at a 180° 
contact angle. Using apriori energy estimates, we derive conditions on variable speed that 
guarantee that a sufficiently smooth solution of the linear advection-diffusion equation blows 
up in a finite time. Using the class of self-similar solutions to the linear advection-diffusion 
equation, we find the blow-up rate of singularity formation. This blow-up rate does not agree 
with previous numerical simulations of the model problem. 

1 Introduction 

Contact lines are denned by the intersection of the rigid and free boundaries of the flow. Flows 
with the contact line at a 180° contact angle were discussed in [21 [6], where corresponding solutions 
of the Navier-Stokes equations were shown to have no physical meanings. Recently, a different 
approach based on the lubrication approximation and thin film equations was developed by 
Benilov & Vynnycky pQ. 

As a particularly simple model for the flow shown on Figure [H the authors of [1] derived the 
nonlinear advection-diffusion equation for the free boundary h(x,t) of the flow: 



dh d 
dt dx 



h 3 ( ,d 3 h dh\ rr , ' 



0, x > 0, t > 0, 



(1.1) 



where a is a numerical constant. The contact line is fixed at x = in the reference frame moving 
with the velocity —V(t) and is defined by the boundary conditions h\ x= o = 1 and h x \ x= o = 0. 
The flux conservation gives the boundary condition for h xxx \ x= o = — ^j. For convenience, we 
can fix a 3 = 3. Existence of weak solutions of the thin-film equation (|l.ip for constant V(t) and 
Neumann boundary conditions on a finite interval was recently constructed by Chugunova et al. 

mi]. 

Using further asymptotic reductions with 



h-l = 0(\V\- 1 ), x = 0(\V\~ 1 ' 3 ), t = 0(\V\-^ 3 ), as \V\-kx>, 



(1.2) 
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Figure 1: Schematic picture of the flow between rigid boundaries. 

the authors of [T| reduced the nonlinear equation (jl.ip with a 3 = 3 to the linear advection- 
diffusion equation: 

dh d A h T ,, . dh . 

subject to the boundary conditions 

h\ x= Q = 1, h x \ x= Q = 0, /ii ra |i=o = — -, t>0, (1.4) 

We assume that h, h x , h xx — > as x — > oo: in fact, any constant value of h at infinity is allowed 
thanks to the invariance of the linear advection-diffusion equation (|1.3|) with respect to the shift 
and scaling transformations. Indeed, if h(x,t) solves the boundary- value problem (|1.3|) — (|1.4|) . 
then 

h(x,t)=H+(l-H)h(£,T), £ ' ' 



(\-H)y^ (i-^) 4 /3' 

with constant H < 1 solves the same advection-diffusion equation (|1.3|) with the same boundary 
conditions (jl.4p but for the variable speed V(t) = j^jj and with the asymptotic value h(x, t) — > H 
as x — > oo. 

With three boundary conditions at x = and the decay conditions as x — > oo, the initial- 
value problem for equation (jl.3p is over-determined and the third (over-determining) boundary 
condition at x = is used to find the dependence of V on t. Local existence of solutions to the 
boundary- value problem (jl.3jl - (|1.4jl was proved in our work [7] using Laplace transform in x and 
the fractional power series expansion in t. 

We shall consider the time evolution of the boundary-value problem (|1.3j )- (|1.4j) starting with 
the initial data h\t=o = ho(x) for a suitable function ho. In particular, we assume that the profile 
ho(x) decays monotonically to zero as x — > oo and that is a non-degenerate maximum of ho 
such that ho(0) = 1, h' (0) = 0, and h' \0) < 0. If the solution h(x,t) losses monotonicity in x 
during the dynamical evolution, for instance, due to the value of 

f3(t) := h xx (0,t) (1.5) 



2 



crossing from the negative side, then we say that the flow becomes non-physical for further 
times and the model breaks. Simultaneously, this may mean that the velocity V(t) blows up, 
as it is defined for sufficiently strong solutions of the advection-diffusion equation (jl.3p by the 
pointwise equation: 

h xxxxx {0,t) = V{t)P(t), (1.6) 

which follows by differentiation of (|1.3j) in x and setting x — > 0. 

The main claim of pQ based on numerical computations of the reduced equation (jl.3p as well 
as more complicated thin- film equations is that for any suitable ho, there is a finite positive time 
to such that V(t) — > — oo and (3{t) — > 0~ as 1 1 to- Moreover, it is claimed that V(t) behaves near 
the blowup time as the logarithmic function of t, e.g. 

f(t)~Cilog(to-t) + C 2j as ttto, (1.7) 

where C\, C2 are positive constants. 

The goal of this paper is to inspect possible blow-up rates of the singularity formation in the 
boundary- value problem (|1.3p - (|1.4p . First, we use apriori energy estimates to show that V(t) 
cannot remain positive for all times for smooth solutions of the boundary-value problem (|1.3p ~ 
(jl.4p . This result implies simultaneously two things: if V(t) remains positive, the smooth solution 
blows up in a finite time, and if a smooth solution exists for all times, then V(t) either oscillates 
or become negative. Similarly, we also show that f3(t) and V{t)f3 2 (t) cannot remain negative for 
all times in the same sense: if f3(t) and V{t)j3 2 {t) remain negative, the smooth solution blows 
up in a finite time, and if a smooth solution exists for all times, then f3(t) and V(t)(3 2 (t) either 
oscillate or become positive. Combination of both results shows that the only way a smooth 
solution can exist for all times is if the variable speed V(t) oscillates from positive to negative 
values back and forth. 

Second, we study the class of self-similar solutions based on the scaling transformation (|1.2p . 
The class of self-similar solutions is defined by the linear advection-diffusion equation (jl.3p . 
the decay condition at infinity, and the first two boundary conditions (II. 4p . The third boundary 
condition h xxx \ x= o = — \ is not satisfied for the self-similar solutions and we replace this boundary 
condition with new condition h xxx \ x= o = JoV(t) for a fixed 70 < 0. We show that the solution 
blows up in a finite time for positive V{t) and positive /3(t), which agrees with the scaling 
transformation (ll.2p but does not correspond to the physical requirements of the flow on Figure 

m 

Finally, we study how f3(t) may vanish and V(t) may diverge in a finite time by using the 
pointwise equation (|1.6p and its derivative. We find yet another rate of singularity formulations, 
which is different from the rates based on the scaling transformation (11.2p and on the numeri- 
cally claimed result (jl.7p . Therefore, further studies of the boundary- value problem (jl.3p - (jl.4p 
including more precise numerical studies are required. These studies will be reported elsewhere. 

The remainder of this paper is organized as follows. Section 2 gives apriori energy estimates 
for the boundary-value problem (jl.3p - (jl.4p . Section 3 describes self-similar solutions describing 
blow-up rate of the singularity formulations. Section 4 reports analysis following from pointwise 
equations. 
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2 Apriori energy estimates 



Let us consider the advection-diffusion equation (jl.3p subject to the boundary conditions (|1.4|) 
and the decay condition h, h x , h xx — > as x — > oo. We assume existence of a smooth solution to 
this boundary-value problem and show that V(t) cannot remain positive for all times. 

Theorem 1 Solutions of the boundary-value problem 1 1.3\) - [L4\ ) do not exist in class h G 
_,L 2 (IR + )) nL 2 (R + ,H 2 {R + )) ifV(t) > V > -1 for all t > t > 0. 



Proof. From the advection-diffusion equation (|1.3p . we have the energy balance: 

d t Q/i 2 ) + d x (hh xxx - h x h xx - ^V{t)h 2 ^j + (h xx ) 2 = 0. 

Integrating this equation in x on (0, oo) and using the boundary conditions (jl.4j) and the decay 
conditions as x — > oo, we obtain apriori energy estimates: 



-i 

dt 1 



|2 



) + 2\\h xx (;t)\\h iR+) = -(l + V(t)). (2.1) 



If we have a solution in class h G C(M+, L 2 (M + )) n L 2 (M + , i7 2 (R + )), then integrating the apriori 
energy estimate (|2.ip . we obtain 



\\H;t)\\h { R +) + 2f o ll^(-^)lli 2 (K + )^=ll^o||!2 (R+) -£(l + n^))^- (2-2) 

Since the left-hand-side is strictly positive, the assertion of the theorem is proved. □ 
Next, we rewrite the advection-diffusion equation (I1.3P for the variable u = h x in the form 

ut + u xxxx = V(t)u x , x > 0, t > 0, (2.3) 

subject to the boundary conditions at the contact line 

u| x= o = 0, u^l^o = — -, U xxx |a; = o = 0, t > 0, (2.4) 

where the boundary conditions u xxx \ x= o = h xxxx \ x= Q = follows from the boundary conditions 
h\ x= Q = 1 and h x \ x =o = as well as the advection-diffusion equation (|1.3p as x — > 0. Denote 
(3(t) = h xx \ x= o = u x \ x= o and recall that /3(0) < initially. Again, we assume existence of a 
smooth solution to the boundary- value problem (|2,3p - (|2.4p and show that j3(t) and V(t)f3 2 (t) 
cannot remain negative for all times. 



Theorem 2 Solutions of the boundary-value problem \2.3\) - (K4^ do not exist in class u G 
C{R+,L 2 {R + ))nL 2 (R + ,H 2 {R + )) if 0(t) < /3 < for all t > t Q > 0. 
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Proof. From the advection-diffusion equation (|2.3p . we have the energy balance: 
dt + ®* (^ uu xxx ~ UxUxx ~ ^V(t)u 2 ^j + (u xx ) 2 = 0. 

Integrating this equation in x on (0, oo) and using the boundary conditions (|2.4p and the decay 
conditions as x — > oo, we obtain apriori energy estimates: 

d 11 ' — 12 , oil, r_ +MI2 



dt 



<Mh(R + )+ 2 \\ u U-,t)\\h(R + ) = £(*)■ ( 2 - 5 ) 



If we have a solution in class u G C(R+, L 2 (R+)) R L 2 (R+, if 2 (R+)), then integrating the apriori 
energy estimate (12. 5h . we obtain 



ft 

IH-> t )llL2(R + ) + 2 ^ W u xx(-,T)\\ 2 L 2 {R+) dT = \\u \\l 2{R+) + p(r)dr. (2.6) 
Since the left-hand-side is strictly positive, the assertion of the theorem follows. □ 



Theorem 3 Solutions of the boundary-value problem 112. ,3\) - (K4\ ) do not exist in class u G 
-,H 1 (R + ))nL 2 {R + ,H 3 (R + )) ifV{t) < V < and \p(t)\ > /3 > for all t > t > 0. 



Proof. Multiplying the advection-diffusion equation (|2 . 3[) by u xx , integrating this equation 
in x on (0, oo), and using the boundary conditions (|2.4|) and the decay conditions as x — > oo, we 
obtain apriori energy estimates: 

^IM-,t)|| 2 2(R+) + 2||u :raa; (-,t)|| 2 2{R+) = V{t)f3 2 (t). (2.7) 

If we have a solution in class u G C(R+, -ff 1 (M+)) n L 2 (R + , iJ 3 (R + )), then integrating the apriori 
energy estimate (|2.7f) . we obtain 

t ft 

2 J ll„./ ||2 , / Trf„\o2i 



\\M-,t)\\h(R + ) + 2 J o hxxx(;r)\\i 2{1&+) dT= K||£ a(R+) + y F(t)/3 2 (t)cZt. (2.8) 
Since the left-hand-side is strictly positive, the assertion of the theorem follows. □ 

3 Self-similar solutions for singularity formations 

Let us consider the class of self-similar solutions to the linear advection-diffusion equation (jl .3j> : 



where to is an arbitrary positive parameter for a finite blowup time, Vq is an arbitrary parameter 
for the initial velocity, and /(£) is a solution of the differential equation: 



d 4 f 1 r/f 

|/ + _ (C _4t Vb)| = 0, e>0. CJ.2) 
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We are looking at a solution of the boundary- value problem associated with the first two boundary 
conditions at the contact line: 

/(0) = 1, /'(0) = 0, (3.3) 

and the decay condition /(£), /'(£) — > as £ — > oo. Note that the third condition at the contact 
line h xxx \ x= Q = — \ is not satisfied by the self-similar solution (|3.ip . The revised third boundary 
condition is given by 

f"'(0) 

h xxx \ x=0 = (tQ _ t) 3 /4 = ioV(t) , (3.4) 

where 70 is constant such that f"'(0) = io^oTo- Also note that the class of self-similar solutions 
(|3.ip is compatible with the asymptotic scaling (jl.2p used in the derivation of the linear advection- 
diffusion equation (|1.3p . 
Setting 

g(z) = f'(0, z = £-4t V , 
we reduce the boundary-value problem (|3.2p — (|3.3|) to the following system: 

Ag"'(z) + zg{z) = 0, z > z , 

g{z ) = 0, (3.5) 
g(z), g'(z) — > as z — » 00, 

where zq = — 4toVo- A suitable solution of this boundary-value problem is constructed in the 
following theorem. 

Theorem 4 There exists a unique (up to scalar multiplication) positive solution of the boundary- 
value problem 113. 5|) on (z ,oo) with zq < 0. 

Proof. As z — > 00, there are three fundamental solutions of the linear equation 

Ag"'(z)+zg(z) = 0. (3.6) 

One solution decays to monotonically as z — > 00 and the other two solutions oscillate and 
diverge as z — > 00. Therefore, the space of solutions of the boundary-value problem (|3.5j) is 
spanned by a particular solution (denoted by G) decaying to at infinity. To define G uniquely, 
we construct a decaying solution of the differential equation (|3.6p asymptotically by using the 
WKB analysis 0: 



G(z) = exp 



3z 4 / 3 ' 
" 2 8 / 3 



+ 



1/3 V z 2/3 



as z — > 00, (3-7) 



where corrections terms can be identified in terms of power series in inverse powers ofz 1 / 3 . The 
solution G of the linear equation (|3.6p can be extended globally for all z G R. To satisfy the 
boundary condition at zq, it remains to show that there is zq £ R such that G(zq) = 0. 

It is clear that zq € R exists. Indeed, if zq does not exist, then G{z) remains positive for 
all z S R, which is only possible if G(z) decays to monotonically as z — > —00 (the other two 
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solutions again oscillate and diverge as z — > — oo). However, then G G H 2 (M.) is a global solution 
of the differential equation (|3.6[) for all z£l. Multiplying this equation by G' and integrating 
by parts, we obtain the contradiction 

4 f {G"fdz + - [ G 2 dz = 0, (3.8) 

which proves that no G G H 2 (M.) may exist. Furthermore, zq < because G(z) is monotonically 
decaying for all z > 0. To see this, we use the fact that the differential equation (|3.6p is invariant 
under the transformation z \— > —z, so that G{z) := G(—z) is another solution of (13. 6p . The func- 
tion G{z) increases monotonically for large negative z. Since 4G'"(z) = —zG(z) > for all z < 0, 
G(z) remains monotonically increasing for all z < and hence G(z) decreases monotonically for 
all z > 0. Therefore, zq < 0, that is, Vo > (if to > 0). The value zo is uniquely determined as 
the largest negative zero of the positive function G(z). □ 
Figure [2] shows numerical approximation of the solution G{z) satisfying the boundary-value 
problem ()3.2p . The numerical approximation is obtained with the standard Heun method. 




Figure 2: Numerical approximation of the solution G(z) of the boundary-value problem (|3.2j) . 

A general solution of the boundary- value problem (|3.5h is given by g(z) = CG(z). To deter- 
mine the constant C, we use the modified boundary condition (|3.4p . If h xxx \t=o < as inherited 
from the third boundary condition in (jl.4p . we should have g"(0) = f"'(0) < or CG"(zq) < 0. 
By the continuity arguments, we have G"(zq) < and therefore, C > 0. Indeed, G(z) := G(—z) 
is monotonically increasing function for all z < with AG"'{z) = —zG(z) > 0. When z > 0, 
G"'{z) < as long as G(z) > 0, so that there is zq G (0,2^), such that G"{z) < for all 
z G (— Jo — ^o]) or equivalently, G"{z) < for all z G [zo, zo)- Therefore, G"(zo) < as on Figure 

El 

By the same argument, there is zq G (zo,zo\ such that G'(z) < for all z G (— zq,— zq[, or 
equivalently, G'(z) > for all z G [zo^o)- Therefore, G'(zq) > as on Figure El which implies 
that 

/"(0) CG'(zo) r \ 
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We conclude that the class of self-similar solutions (|3.ip does not represent the relevant 
dynamics of the boundary-value problem (|1.3j )- (|1.4p in the context of the physical requirements 
of the flow on Figure [T] because /3(t) is supposed to be negative at least for initial values of t > 0. 

4 Point wise equations 

We give here additional estimates of how the solution of the boundary-value problem (11.3j) - (11.4p 
may blow up in a finite time, based on the pointwise equation (| 1 .6 j) and its derivative. We look at 
the boundary- value problem (|2.3p - (|2.4p and assume existence of a sufficiently smooth solution. 
By taking the limit x — > 0, we recover the pointwise equation (|1.6p rewritten in new variables as 

V(t)f3(t), t > 0, (4.1) 

where f3(t) = u x \ x= q. By taking a derivative of the linear advection-diffusion equation (|2.3[) in x 
and the limit x — > 0, we obtain another pointwise equation: 

^Tf + u xxxxx \x=o = ~^V(t), t > 0. (4.2) 

The system of equations (|4,ip and (|4.2p can be rewritten in the partially closed form: 

dfi u xxxx \ x= q . , . 

~dt ~ WW ~ Uxxxxxlx=0 ' - u - 1 6) 

Let us now assume that there is to > such that 

Pit) — s- 0, u xxxx \ x= o — > 04, u xxxxx \ x= q — > as, as t f to, (4.4) 

where ^ and |as| < oo. Then, asymptotic analysis of the differential equation (|4.3p shows 
that 

P 2 {t) = a A {t Q -t) + 0{to-tf/\ V(t) = J-^- t +G(l), as t\t , (4.5) 



under the constraint that a A > 0. The asymptotic rate (I4.5h is different both from the scaling 
transformation (|1.2p and the numerically claimed result fjl .7|) . In the context of the numerical 
result (jl.7p . this pointwise analysis may imply that either = or 05 — >■ 00 in the assumption 

We conclude that three different rates of the singularity formations claimed in (II. 7p and 
obtained in (|3.ip and (j4.5j) indicate complexity of dynamics of the boundary- value problem (|1.3j) 
(ll.4p or its equivalent version (I2.3p ~ (l2.4p . Further studies of dynamical evolution of contact lines 
within this reduced problem are needed, including more precise numerical simulations. 

Acknowledgement: The authors are thankful to E.S. Benilov and R. Taranets for useful 
discussions and for sharing their unpublished results at an early stage of his research. 
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